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ABSTRACT 

Context. The elliptical galaxy NGC 3923 is known to be surrounded by a number of stellar shells, probable remnants of an accreted 
galaxy. Despite its uniqueness, the deepest images of its outskirts come from the 1980s. On the basis of the modified Newtonian 
dynamics (MOND), it has recently been predicted that a new shell lies in this region. 

Aims. We obtain the deepest image ever of the galaxy, map the tidal features in it, and search for the predicted shell. 

Methods. The image of the galaxy was taken by the MegaCam camera at the Canada-France-Flawaii Telescope in the g' band. It 
reached the surface-brightness limit of 29 mag arcsec“^. In addition, we reanalyzed an archival HST image of the galaxy. 

Results. We detected up to 42 shells in NGC 3923. This is by far the highest number among all shell galaxies. We present the 
description of the shells and other tidal features in the galaxy. A probable progenitor of some of these features was discovered. The 
shell system likely originates from two or more progenitors. The predicted shell was not detected, but the new image revealed that the 
prediction was based on incorrect assumptions and poor data. 

Key words. Galaxies: peculiar - Galaxies: elliptical and lenticular, cD - Galaxies: individual: NGC 3923 - Galaxies: photometry - 
Techniques: image processing - Gravitation 


1. Introduction 

Stellar shells in elliptical galaxies are likely remnants of g alactic 
interactions (see |Dupraz & Combes|1986[|Ebrova|2()T3[ or |Bflek| 
et al. [2015] for a review). The most accepted formation scenario 
is the accretion of a small satellite along an almost radial trajec¬ 
tory ( |Quinn|19'8^ where the stars forming the shells are released 
from the progenitor when it passes through the pericenter. These 
stars form the shells when they reach their apocenters. For the 
Quinn model [Pupraz & Combes] ( p^87) ) showed that the progen¬ 
itor has to make a few oscillations before it dissolves to explain 
the high radial range (i.e., the ratio of the radii of the outermost 
and the innermost shell) observed in many shell galaxies. The 
accreted galaxy loses orbital energy by dynamical friction. The 
outermost shells are created first and the innermost shells later 
when the progenitor loses enough energy. Shells that were cre¬ 
ated together during a particular pericentric passages are called 
a generation (see |Bilek et al.|2013| for details). 

( |1984| ) pointed out the connection of the shell 


Quinn 


radii with the gravitational potential of the host galaxy. |Bflek| 


et al.| ( [2014| ) (hereafter B14) predicted the existence of an 


as yet unobserved shell in the elliptical galaxy NGC 3923 at 
-2000"southwest (around 200 kpc) from the galaxy’s center. 
The prediction wa s made supposin g the modified Newtonian dy¬ 
namics (MOND) ( |Milgrom|1983| ), a modification of Newtonian 
dynamics for low accelerations that eliminates the need for dark 
matter. This prediction was based on the shell radii tabulated by 


|Prieur| ( |1988| ) and [Sikkema et al.| ( |2007| ). If the shell were to 
be discovered, it would be the first discovery of an object pre¬ 
dicted by MOND, a test of MOND to exceptionally low gravita¬ 
tional acceleration for an elliptical galaxy, and the discovery of 
the biggest shell ever. 

Here we present the new deep observation of the galaxy we 
obtained to search for the predicted shell. The observations of 
the outer parts of the galaxy (above a radius of 140") supersede 
the best existing images so far coming from the 1980s. More¬ 
over, we reanalyze the HST image obtained by [Sikkema et al.| 
( |2006 ) which we use to look for the shells in the inner parts (un¬ 
der 140") of the galaxy. We present a new table of shell radii 
and describe the other tidal features in the galaxy that can serve 
to constrain the gravitational potential (e.g., |Bflek et al.|[2013| ). 
The predicted shell is not detected, but we find that the analy¬ 
sis of B14 was based on incorrect assumptions and poor data. In 
a forthcoming paper, we will present a quantitative photometric 
analysis. 

The distance of the galaxy from Earth is between 19.9 to 
24 Mpc according to the direct measurements in the NAS A/IPAC 
Extragalactic Database. This corresponds to a linear scale of 
around 1 kpc per 10". 

The paper is organized as follows. In Sect. we describe 
observations, the data reduction, and the methods used to en¬ 
hance fine structures present in the galaxy. These structures are 
described in Sect. Implications of the observation are dis¬ 
cussed and concluded in Sect.|^ Appendix |A| contains detailed 


Article number, page 1 of 14 



















































A&A proofs: manuscript no. N3923_Mega 


images of all shells and other features found in the galaxy. In 
Appendix]^ we describe our minimum masking method for the 
shell detection in detail. 


2. Observations and data reduction 

Table 1: Mosaic tile centers 


Tile 

R.A. 

Dec. 

1 

11:49:30 

-29:11:00 

2 

11:52:30 

-28:25:00 

3 

11:49:30 

-28:25:00 

4 

11:52:30 

-29:11:00 


2.1. New observations by MegaCam 

Our new observations of NGC 3923 were carried out using the 
MegaCam imager at the 3.6 m Canada-France Hawaii Telescope 
(CFHT) within the program 14B012. They took place in January 
and February 2015. The galaxy was covered by a four-tile mo¬ 
saic (see Tab. with the resulting field of view of around 2x2° 
centered on the galaxy. Each tile is a stack of 21 subexposures 
with the exposure time of 195 s in the g' band. The total exposure 
time was 273 min. The pointing offsets between the subframes 
were 7'. These large offsets were required to subtract the para¬ 
sitic light scattered within the instrument. It was necessary to fit 
the imaging of each tile to about 1 hour to avoid substantial vari¬ 
ation of the atmospheric illumination. [Due et al.|p015] ) explain 
in detail the observation strategy needed to capture low-surface- 
brightness features with MegaCam. This camera has the angular 
resolution of 0.185'7pixel. 

The raw data were proces sed by the Elixir pipeline at CFHT 
( [Magnier & Cuillandre|2004| ) to deliver fully detrended images 
(standard CCD processing plus uniform zero-point illumination 
correction, and absolute photometric and astrometric calibra¬ 
tion). Those data were then fed into the low-surface-brightness 
branch of E lixir (Elixir-LSB, |Cuillandre & Ferrarese|20TT| |Duc| 
et al.||2015| ) in order to subtract the large-scale component due 


to parasitic light scattered within the instrument. Elixir-LSB re¬ 
stores the true sky background to a flatness level lower than 0.2% 
of the absolute sky level (min-max), that is more than 7 magni¬ 
tudes fainter, allowing the detection of extended features as faint 
as 29 mag arcsec"^ against a flat background. This limit (stan¬ 
dard for all Elixir-LSB data; see [Ferrarese et al.|[MT^ is es¬ 
tablished through the visual identification of the faintest parts 
of extended features in the image such as shells, tidal streams, 
or galactic cirri. All frames were stacked by MegaPipe (|Cwyn| 

|20fe| ). 


2.2. Older observations by HST 

We studied our MegaCam image together with an image taken 
by the Hubble space telescope (HST) capturing t he center of the 
galaxy. The HST observation was presented by |Sikkema et al.| 
( 2006| ) and |Sikkema et ^.| ( [2007| ). We obtained the data from the 
Hubble Legacy Archiv^We worked with a stack of the images 
taken by the ACS/WFC camera in the filters F814W (978 s) and 
F606W (1140 s). This made the exposure time of 35 min in total. 
The angular resolution of this image was 0.05"/pixel. 

^ http://hla.stsci.edu/ 
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2.3. Shell detection methods 


Shells are typically very faint, low-contrast features (the exact 
values of the surface brightness of the shells will be presented in 
a forthcoming paper). Their visibility benefits from some image 
processing for contrast improvement. We looked for shells via 
a visual inspection of images processed in one of the following 
ways: 

No processing - Most of the reported shells are visible in the 
original FITS images without further processing, although many 
of them with difficulty. 

Model subtraction - We fit th e image of the g alaxy by 
a smooth analytic profile in Calfit ( |Peng et S^|20Q2 ). We use 
a Sersic profile or the sum of two Sersic profiles, both combined 
with sky. The brightest stars, shells, tidal features, and galactic 
cirri are masked out before the fitting. The resulting image is 
the difference between the original image and the model. The 
method is used in Figs. |A.6| and |A.9f[A.12| We use a special fit 
for every image. 

Minimum masking - A median filter is first applied to the 
original image to smooth the data on small angular scales. Next, 
an erosion (minimum) filter is applied. The result is the differ¬ 
ence between the original image and the filtered image. Our ero¬ 
sion filter is circular and the median filter is square. The side 
of the square is the same as the radius of the circle. As far as 
we know, minimum masking has never been described before 
(further details in Appendix This processing does not af¬ 
fect the position of shell edges. Shell edges transform into bright 
filaments on the dark background (in the positive image). Fig¬ 
ures |A!2J|A^ and |A!7}|A^ have been processed by this method. 
In the figures, the term “size” refers to the radius of the erosion 
filter. We tried several filter sizes. This helped us to further con¬ 
firm the existence of the shell, see Fig. |B.5[ 

Unsharp masking - The original image is convolved with 
a Gaussian. The result at every pixel is twice the original inten¬ 
sity minus the intensity in the convolved image. We use the stan¬ 
dard deviations of the Gaussian of 44.4, 22.2, and 11.1" for our 
MegaCam image. This method proved to be less effective than 
minimum masking and model subtraction since unsharp masking 
produces dark rings around bright stars that disturb the images 
and complicate the detection of the shells. This method helped 
us only to confirm the existence of some shells. 

Before including a shell in our list (Tab.[^, it had to be visi¬ 
ble in images processed by two different methods at least. 


We detected more shells in the HST image than Sikkema 
|et al.| ( |2Q07| ) because they used only the model subtraction 
method. Some shells become more pronounced by the minimum 
masking method. On the other hand, while the model subtraction 
clearly reveals, for example, shell S40 (Figs. |A.6] and |A.ll]), this 
shell is barely visible after the minimum masking (Fig. |A.7| ). 


3. Results 

3.1. Shells 

The shells or shell candidates we found are described in Tab. 
where we present the radius of the shell in arcseconds, its promi¬ 
nence degree, an image processing method allowing its reliable 
detection, the data the shell was found in, the enumeration of the 
images showing shell, and notes. The shells were detected by eye 
and thus the prominence rating is subjective. We state only the 
image processing method allowing a reliable detection with the 
minimal effort (order of difficulty: no processing, unsharp mask¬ 
ing, minimum masking, model subtraction). Seeing a shell in an 




















































M. Bilek et al.: Deep imaging of the shell elliptical galaxy NGC 3923 with MegaCam 



Fig. 1: Shell appearance differs depending on the line of sight. In the three-dimensional space, a shell is a part of a spherical surface 
centered on the galaxy. The observer in direction A sees the shell as a sharp-edged structure. The observer viewing the shell from 
direction B does not see any part of the shell tangentially, and so the shell appears diffuse. The two images on the right show an 
illustrative simulation of a radial minor merger between two spherical galaxies. Only the stars from the smaller galaxy are shown. 
The middle panel shows the view perpendicular to the line of collision. The view in the right panel is inclined by 35° to the line of 
collision. 


image processed by one method, we could note it in almost all 
the images processed in other ways. The shell system is approx¬ 
imately axially symmetric. The symmetry axis coincides with 
the major photometric axis of the galaxy. Their position angle is 
close to 225°(the southwest-northeast direction). 

The shells that have the radii with the plus signs in Tab. lie 
on the northeastern half of the galaxy and those with minus signs 
on the southwestern half. We measured the radius of every shell 
along the major photometric axis of the galaxy. We measured the 
radius from the brightest pixel in the center of the galaxy to the 
outer cut-off of the shell. The assumed galaxy center had the 
FK5 coordinates of R.A. = 11:51:01.7 and Dec. = -28:48:21.55. 

In a forthcoming paper, the table will be completed by the 
shell equatorial coordinates, surface brightness estimates, and 
the shell azimuthal widths. This will be interesting for compari¬ 
son with simulations. 

In the center, we looked for the shells independently in the 
MegaCam and the HST image. The shells were more visible in 
the HST image because of its higher angular resolution and the 
absence of seeing. As can be noted from Tab.[^ we missed sev¬ 
eral shells in the MegaCam image. All of them are seen as very 
diffuse in the HST image. On the contrary, we detected only one 
shell candidate in the MegaCam image that we missed in the 
HST image (radius of -87"). It would have been rated at the low¬ 
est prominence degree. In the HST image this feature looks like 
a short filament rather than a shell and so we do not include it 
into our table. For the shells captured in the HST image, the table 
states the radii measured from this image. 

All the shells we detected are shown in Figs. |A.l A.7 They 
are marked by S^. 

3.2. Other tidal features 


Figures |A.8f[A.12| were processed to show the other tidal fea¬ 
tures in the galaxy. 

Shell S3 has a much narrower azimuthal extent than the other 
shells in the galaxy. We can see that it lies at the end of a stream 


likely the progenitor of the shell S3 and possibly of some oth¬ 
ers. The projected distance of this galaxy from the center of 
NGC 3923 is 474" (around 47 kpc supposing that it lies close 
to NGC 3923). It is the elliptical g alaxy PGC 3097920. A ccord- 
ing to the HyperLed^ database ( Makarov et al. 2014), it has 
a heliocentric radial velocity of 1663+39 km s“\ while that of 
NGC 3923 is 1550+15 km s“\ and so they can be in physical 
proximity. 

A hook-like structure (H in Figs. A.l[ |A.8||A.9| and |A.10| ) 
lies northwest of the galaxy center. First, it goes to the north from 
the galaxy center, then it bends to the west and disappears point¬ 
ing to the south. Similar structures can be produced by an accre- 
tion of a disk galaxy as we can see in the simulations depicted in 
|Dupraz & Combes] ( |1986| ). On the other hand, this feature could 
also be a galactic cirrus. The western end of the hook is located 
near a triangular object (Cl in the hgures) with a hlamentary 
structure typical for cirri. The hook itself appears smooth. 

In the neighborhood of the hook and the shell S5, there is 
a fan feature (F in Figs. |A.l[|A.9| and |A.10| ). It is reminiscent of 
a shell, but its edge is not tangential to the direction toward the 
galaxy center and it has a lower azimuthal extent than most of 
the shells in the galaxy. 

Another prominent tidal feature is the stream (labeled ST in 
Figs. [ATljA3] and |A:8HAT0l ) south of the core of NGC 3923. It 
points almost precisely in the north-south direction. It has no 
internal structure. It does not seem to be an extension of the 
hook (Fig. \A.^ . It disappears near the edge of the shell S4 on 
the south. At hrst sight, it appears to be a reflection caused by 
the star below it. However, if it were an artifact, it would have 
to move with respect to the stars during our large offsets, but 
when we checked the individual subexposures, it was not the 
case. Moreover, it is too fuzzy to be a common artifact. 

As has been done in previous observations (e.g., [Sikkema 
|et al.||2007] |Pence|p^^ , we note several dust featur es nea r 
the center of the galaxy in our MegaCam image (Fig. |A.ll ) 
Pence] ( |1986|) r eported that the most prominent cloud (labeled 


extending approximately from the core of NGC 3923 (Figs. A.9 
|A.10[ and |A.12| ). A compact galaxy lies on the symmetry axis 


DP in Fig. |A.4| ) contains Ha emission with the same redshift as 
NGC 3923. This proves that the dust lies in the galaxy and not in 
the foreground. 


of this stream. The stream gets thinner close to the galaxy. It is ^ http://leda.univ-lyonl.fr/ 
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3.3. Galactic cirri 

We can see many galactic cirri in our MegaCam image (labeled 
C^). They are dust clouds in our own Galaxy. They typically 
show a fine filamentary structure. As is well known, galactic cirri 
can easily be mistaken for tidal features in the external galaxies 
even if we have observations in several filters (e.g., [Due et aL| 
|2015| ). The angularly large cirri can be identified by comparing 
the optical images with the radiation maps in 857 GHz obtained 
by the Planck spacecraft ( |Duc et al.|2015| ). There is a structure 
in the near-infrared map which could correspond to the stream. 
We did not identify such a counterpart for any of the other tidal 
features we report, but the possibility of confusion is not fully 
excluded. 


4. Discussion and conciusions 


We presented our new deep image of the shell elliptical galaxy 
NGC 3923. It was taken in the g' band and reaches the surface- 
brightness limit of 29magarcsec“^. In addition, we reanalyzed 
an archival HST image of the galaxy. The motivation for our new 
observation was to look for the shell predicted by B14 to make 
a test of the MOND theory of modified dynamics. The predicted 
shell was not found. 

We detected 42 shells or shell candidates in total. Even 
though a few candidates are highly doubtful, the number is sur¬ 
prisingly higher than previously reported. The compilation by 
B14 contains only 27 shells. Two of the shells listed th ere al 
legedly lying at the radii -r630"and -520" do not exist. IPrieur 


( |1988| ) possibly confused the halves of the galaxy. The shell at 
-r630" is probably S4 with the radius of -630" and the shell at 
-520" is perhaps S5 with the radius of +490". This means that 
the prediction of B14 was based on incorrect and incomplete 
data. Next, this teaches us that many shells (~17) can easily es¬ 
cape observations, not a maximum of only six as assumed by 
B14. 

A probable progenitor of at least one of the shells was found. 
It is the galaxy PGC 3097920. Since the ratio of the projected 
radii of the innermost prominent shell S41 and this satellite is 26, 
the shell system was likely created by two or more progenitors. 
This follows from the findings from simulations summarized in 
|Dupraz & Combes] ( |1987| ). It was found that the shell systems 
with a radial range higher than around 6 had to be made of mul¬ 
tiple generations of shells (see Sect. in the framework of the 
minor merger model. The largest shells form first, while the in¬ 
ner shells form later when the secondary loses enough orbital 
energy by dynamical friction. But even if the observed progeni¬ 
tor is now in the apocenter, it has enough energy to create shells 
as large as some of the largest shells in NGC 3923 in its next 
pericentric passage. Thus, it seems that the shells had to be cre¬ 
ated by more than one progenitor. This is supported by the fact 
that the number of shells in NGC 3923 is the highest of all galax¬ 
ies. The works of |Bilek et al.| ( |2013] ) and B14 assumed that the 
whole system was created from one progenitor. 

For these reasons, our observation has no implication for the 
validity of MOND at this moment. The consistency of the shell 
distribution with MOND [Bflek et ( |2013| ) will have to be re¬ 
analyzed. 

The shell system in NGC 3923 held several records before 
our observation; it contained the largest known shell, the high¬ 
est number of shells, and it had the highest radial range. After 
our observation, the galaxy becomes even more exceptional. We 
increased the number of shells and the radial range raised from 
65 to 108. There are several reasons for the substantial increase 
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in the number of shells: 1) the depth of our MegaCam image, 
2) other image processing methods applied to the HST image, 
and 3) more emphasis on the complete detection of shells moti¬ 
vated by the possibility of constraining the gravitational potential 
dBflek et al.|20T3l ). 

We note in the figures and Tab.j^that the well-defined shells 
are mixed among the diffuse and very faint shells. There are sev¬ 
eral reasons for this. For example, the shells can come from dif¬ 
ferent shell generations or from different progenitors. It is also 
possible that the diffuse shells are shells seen from unfavorable 
angles. Indeed, shells are parts of spheres in three-dimensional 
space. If the shell is to be observed as a sharp-edged feature, the 
line of sight must be tangential to the shell surface (Fig. [T]). In 
the opposite case, the shell edge appears diffuse or the shell is 
not observable. 
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Table 2: Detected shell candidates in NGC 3923 


Label 

Major-axis 
radius ["] 

Prominence 

Data 

Detection 

method 

Figures 


Note 

SI 

+ 1170 

2 

MC 

1 

A.l 



largest 

S2 

-952 

4 

MC 

1 

A.l 

A.2 


MP, 907-1012", highly uncertain, possibly merged 
star halos 

S3 

-846 

1 

MC 

1 

A.l 

A.2 

A.8 

narrow, connected with the “progenitor” 

S4 

-630 

3 

MC 

1 

A.l 

A.2 

A.8 

MP, 640-590", irregular, diffuse, near the end of the 
“stream” 

S5 

+490 

3 

MC 

1 

A.l 

A.3 

A.8 

MP, possibly other tidal feature, small azimuthal 
extent, opposite S8, connected with the “hook”, 
shells at lower radii have surface brightness in¬ 
creased at the azimuth of this shell 

S6 

-430 

4 

MC 

3 

A.2 

A.3 


MP, very diffuse and faint, highly uncertain 

S7 

-1-363 

1 

MC 

1 

A.l 

A.3 

A.8 


S8 

-338 

3 

MC 

3 

A.2 

A.3 

A.8 

MP, not on axis, only in the eastern half, opposite 
S5, shells at lower radii are brighter on the azimuth 
ofS8 

S9 

+331 

3 

MC 

3 

A.3 



copies the edge of S7 but lower azimuthal extent 

SIO 

-280 

1 

MC 

1 

A.3 

A.8 



Sll 

+260 

3 

MC 

3 

A.3 



MP, very diffuse, a smooth step in surface bright¬ 
ness, possibly other tidal feature 

S12 

-237 

4 

MC 

3 

A.3 



MP, sharp-edged, small azimuthal extent, possibly 
other tidal feature 

S13 

+203 

1 

MC 

1 

A.3 

A.8 



S14 

-150 

1 

MC 

1 

A.3 




S15 

+ 146 

2 

MC 

1 

A.3 

A.4 



S16 

+ 130 

1 

MC, HST 

1 

A.4 




S17 

-120 

3 

MC, HST 

1 

A.4 



MP 

S18 

-103 

1 

MC, HST 

1 

A.4 




S19 

+99.5 

3 

MC, HST 

3 

A.4 




S20 

-92.0 

4 

HST 

3 

A.4 



MP, very faint and diffuse, smooth step in surface 
brightness, possibly other tidal feature 

S21 

+90.6 

3 

MC, HST 

3 

A.4 



MP 

S22 

-78.4 

2 

MC, HST 

1 

A.4 

A.5 


sharpness varies with azimuth 

S23 

+72.8 

2 

MC, HST 

1 

A.4 

A.5 

A.6 


S24 

-71.4 

3 

MC, HST 

3 

A.4 

A.5 


MP, not on axis, only in the western half, diffuse, 
possibly other tidal feature 

S25 

-68.6 

1 

MC, HST 

1 

A.4 

A.5 



S26 

+64.8 

3 

MC, HST 

3 

A.4 

A.5 



S27 

+60.0 

2 

MC, HST 

1 

A.4 

A.5 




Notes. Label - Shell designation. Radius - Shell radius in arcseconds. Shells lying on the north-eastern side of the galaxy have the positive sign; 
the shells on the south-western side of the galaxy have the negative sign. Prominence - I - Prominent and sharp-edged, 2 = Prominent and 
diffuse-edged, 3 = Faint but prob ably existent, 4 = Que stionable. Data - The image where the shell was detected: our observation with MegaCam 
(MC) or the archival HST image ( |Sikkema et al.|200^. De tection method - A method allowing a good detection of the shell: 1 = No processing, 
2 = Model subtraction, 3 = Minimum masking (Sect. Appendix [B]). Figures - The figures showing the shell. Notes - MP = Missed by the 
previous observations. 
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Table 2: Continued 


Label Major-axis Prominence 

Data 

Detection Figures 

Note 

radius ["] 


method 



A.4 


A.5 


MP, diffuse, possibly other tidal feature 


S28 

-60.0 

3 

MC, HST 

3 

S29 

-55.9 

1 

MC, HST 

1 

S30 

-r51.0 

1 

MC, HST 

1 

S31 

-49.0 

3 

MC, HST 

3 

S32 

-r45.4 

3 

HST 

3 

S33 

-44.7 

1 

MC, HST 

1 

S34 

-h38.0 

3 

MC, HST 

3 

S35 

-37.5 

3 

MC, HST 

3 

S36 

-r34.8 

1 

MC, HST 

1 

S37 

-30.8 

3 

HST 

2 

S38 

+29.5 

1 

MC, HST 

1 

S39 

-24.2 

4 

MC, HST 

3 

S40 

-rl9.8 

1 

MC, HST 

2 

S41 

-18.1 

1 

MC, HST 

1 

S42 

-10.8 

3 

MC, HST 

3 


A.5 

A.6 



A.5 

A.5 



very prominent 

MP, diffuse, possibly other tidal feature 

A.5 



MP, diffuse, probably the shell at 41.5"in Sikkema 




etal. (2007( 

A.5 

A.6 

, A.7 

very prominent 

A.5 



MP, very faint on axis, well defined on the western 
half 

A.5 

A.7 


diffuse 

A.5 

A.6 

A.7 


A.6 

A.7 


bright smooth step in surface brightness, possibly 
other tidal feature, probably the shell at -28.7"in 
Sikkema et al. (2007|) 

A.5 

A.6 

A.7 


A.6 

A.7 


MP, highly questionable, faint, small azimuthal ex¬ 
tent, possibly edge of a faint dust cloud 

A.6 

A.7 


very prominent after model subraction, abrupt cut¬ 
off in azimuth 

A.6 

A.7 

A.7 


abrupt cutoff in azimuth 

MP, very faint 


Notes. Label - Shell designation. Radius - Shell radius in arcseconds. Shells lying on the northeastern side of the galaxy have a positive sign; 
the shells on the southwestern side of the galaxy have a negative sign. Prominence - I - Prominent and sharp-edged, 2 = Prominent and diffuse- 
edged, 3 = Faint but probab ly existent, 4 = Ques tionable. Data - The image where the shell was detected: our observation with MegaCam (MC) 
or the archival HST image ( |Sikkema et al.|200^. De tection method - A method allowing a good detection of the shell: 1 = No processing, 2 = 
Model subtraction, 3 = Minimum masking (Sect. ]2.3] Appendix]^. Figures - The figures showing the shell. Notes - MP = Missed by the previous 
observations. 


Appendix A: Images 


Here we present images of the individual features we found in 


our MegaC am image and in the archival HST image ( [Sikkema 
et al.||2006| ). In all images, north is to the top and east is to the 


left. The abbreviations used in the images are as follows: 


Sn - Shells. 

A - Artifact caused by reflection of a bright star. It is present 
only on some of the subexposures. 

ST - Stream ( see S ect. |3.2| ). It is present on all subexposures. 

F - Fan (Sect. [ 3 : 2 ]). 

H - Hook (SectjT^. 

DP - Dust patch. 

P - Progenitor of shell S3 and possibly of some others. 

Cn - Probable galactic cirri. They show filamentary structure. 
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M. Bilek et al.: Deep imaging of the shell elliptical galaxy NGC 3923 with MegaCam 



Fig. A.l: MegaCam image. No processing. 




Fig. A.2: MegaCam image. Minimum masking with a filter size 
of 44.4". 


Fig. A.4: HST image. Minimum masking with a filter size of 4". 
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Fig. A.5: HST image. Minimum masking with a filter size of 2". Fig. A.7: HST image. Minimum masking with a filter size of 2". 



Fig. A.6: HST image. Model subtraction. 


Fig. A.8: The area of the hook. MegaCam image. Minimum 
masking with a filter size of 22.2". 
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M. Bilek et al.: Deep imaging of the shell elliptical galaxy NGC 3923 with MegaCam 




Fig. A. 11: Dust in the central parts of NGC 3923. MegaCam im¬ 
age. Model subtraction. 


Fig. A.9: Tidal features around NGC 3923. MegaCam image. 
Model subtraction. 




Fig. A. 12: Progenitor of some tidal features in NGC 3923. Mega¬ 
Cam image. Model subtraction. 


Fig. A. 10: Outskirts of NGC 3923. MegaCam image. Model sub¬ 
traction. 
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Appendix B: Minimum masking 


Here we describe our minimum masking method for detecting 
the shells. The method is based on the operation known as ero¬ 
sion in image processing. The erosion operation is characterized 
by our choice of the neighborhood of every pixel which has a cer¬ 
tain shape and size (rectangle, circle, line, etc.) and is the same 
for all the pixels in the image. At a given pixel, the result of 
erosion is the minimum over its neighborhood in the input im¬ 
age. We used a circular neighborhood. To suppress the effects 
of noise, we smoothed the original image by the median filter 
before applying the erosion filter. For any pixel, the result of the 
median filter is the median over a neighborhood of the pixel. Our 
median filter was square. The side of the square was the same 
as the radius of the circle we used for the minimum filter. This 
proved to be a simple strategy that worked well. In summary, the 
minimum masking method means 1) using a median filter for 
the input image, and 2) applying erosion to the median filtered 
image; the result is the difference between the input image and 
the eroded image. The method is demonstrated in Fig. |B.1| on 
a one-dimensional example (the width of the erosion filter was 
twice the width of the median filter). 


Minimum masking is closely related to the “morphologi¬ 
cal gradient” methods of edge detection in image processing. 
These methods are based on the difference between the origi¬ 
nal image and the eroded or dilated image. The dilatation op¬ 
eration is defined like the erosion, only the minimum is sub¬ 
stituted by the maximum. Moreover, these methods implement 
a noise-smo othing step (e.g., [Gasteratos et al.||1998| or |Yu-qian| 
et al.||200^ . Here the similarity with minimum masking ends. 


The image is then further processed to extract the edges. 


We tested the ability of the minimum masking to detect shells 
in galaxy images. For this purpose, we made a restricted three- 
body simul ation of a shell galaxy formation (see, e.g, |Ebrova| 
|et al.pOl^ for details) by the most widely accepted shell for¬ 
mation mechanism - the phase-wrapping minor merger model 
( Quinn|p~983| ). The two interacting galaxies were modeled as 
Plummer spheres colliding head-on. We stopped the simulation 
when several clearly defined shells developed. The particle po¬ 
sitions were then projected into a plane containing the l ine of 
collision to produce a surface density map (see Fig. |B.6| ). This 
map was used subsequently for our tests. 


In these tests, we were interested only in the rightmost shell 
in Fig. |B.6| to avoid the ambiguity in defining the surface bright¬ 
ness of a shell because the other shells overlap each other in pro¬ 
jection. The chosen shell was made of particles reaching their 
apocenters for the second time since the moment they were torn 
away from their mother galaxy. The radius of this reference shell 
was 135 px (pixels). 


We created two series of artificial images of shell galaxies. 
In the first series, the images were created as a sum of the shell 
image from the simulation and of the image of the body of the 
galaxy modeled by a de Vaucouleurs sphere. The half-light ra¬ 
dius was chosen to be two times the radius of the reference shell. 
Then noise was added in the images. The new intensity of every 
pixel was a choice from the Poisson distribution whose mean 
value was proportional to the original intensity of the pixel. The 
test images in the second series were the same, but the shell im¬ 
age from the simulation was convolved by a Gaussian kernel 
with a standard deviation of 4px. This simulates the observed 
blurriness of shell edges (see the shell surface-brightness radial 
profiles in |Sikkema et al. 12007] ). Each test image is characterized 
by 1) the ratio of the brightest pixel of the reference shell and the 
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Eig. B.l: Demonstration of the minimum masking on a one¬ 
dimensional example. Upper panel: The original data (red) are 
first smoothed by the median filter (green). Then the erosion 
(minimum) filter is applied (blue). Lower panel: The result is 
obtained as the difference between the original and the filtered 
data (the red and blue curve, respectively). The widths of the 
median and erosion filters are indicated in the figure. The unit 
au means arbitrary unit. The original (red curve) is a step-like 
transition between two constant values with added noise. The 
resulting signal (black curve) oscillates around zero outside the 
area indicated by the vertical lines. The width of this area is de¬ 
termined by the width of the transition in the original signal and 
the chosen width of the filter. The original step transforms into 
a bump. 


corresponding underlying pixel of the body of the galaxy (Sh/B), 
and 2) the signal-to-noise ratio (S/N) at the same pixel. 

Einally, we applied minimum masking to the test images. 
The filter size for the first series was 5 px and for the second, 
blurred, series 8 px. These sizes proved to enhance the shell best 
irrespective of the Sh/B or S/N ratios. The results of the first se¬ 
ries are shown in Eig. |B.2| and of the second series in Eig. |B.3[ 
The images show only the area of the shell of interest indicated 
by the red rectangle in Eig. |B.6| 

In the two basic series of test images, the ratio of the ref¬ 
erence shell radius and the characteristic radius of the de Vau¬ 
couleurs sphere was fixed. In Eig. |B.4| we also varied the char¬ 
acteristic radius of the body of the galaxy. We can see that the 
ratio of the shell and the galaxy radius has only a mild influence 
on the shell detectability by the minimum masking method. 

We also compared various methods of shell detection on our 
MegaCam image of NGC3923 to enhance the shell S17 (see 
Eig. |B.5| ). The upper row of tiles shows the image processed by 
the minimum masking method with various filter sizes. The shell 
is clearly detected. Comparing the images processed by the var¬ 
ious filter sizes, we can easily recognize that the shell S17 is not 
an artifact. Tile (d) is processed by a simple linear scaling of the 
original image. Shell S17 is still clearly detectable here. In tile 
(e), a smooth two-component Sersic model of the galaxy body 
was subtracted. The large-scale variations make the detection of 
shell S17 difficult. These variations are caused by the deviations 
of the real galaxy profile from the model. Tile (f) is the image 
processed by unsharp masking. The mask had a Gaussian kernel 
with a standard deviation of 4". This size proved to be the best, 
but even so the shell in question can barely be seen. 
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M. Bilek et al.: Deep imaging of the shell elliptical galaxy NGC 3923 with MegaCam 



Fig. B.6: Simulation of a shell galaxy formation. Only the stars 
from the accreted galaxy are shown. The red rectangle marks the 


region displayed in Figs. B.2-B.4 
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Fig. B.2: Demonstration of the minimum masking method on an artificial image of a shell galaxy. The signal-to-noise ratio S/N 
varies in the horizontal direction. The ratio of the shell intensity to the galaxy body intensity Sh/B varies in the vertical direction. 
The number indicates our rating of prominence. 
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M. Bilek et al.: Deep imaging of the shell elliptical galaxy NGC 3923 with MegaCam 


Signal-to-noise 

100 50 10 5 



Fig. B.3: Demonstration of the minimum masking method on an artificial image of a shell galaxy. The shell component was blurred 
to simulate the observed blurriness of shells. The signal-to-noise ratio S/N varies in the horizontal direction. The ratio of the shell 
intensity to the galaxy body intensity Sh/B varies in the vertical direction. The number indicates our rating of prominence. 


Galaxy-characteristic-radius-to-shell-radius 
1 2 5 10 



Fig. B.4: Demonstration of the minimum masking method on an artificial image of a shell galaxy. The individual images differ by 
the ratio of the shell radius and the half-light radius of the galaxy body. The number indicates our rating of prominence. 
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(d) Linear scaling (e) Model subtraction 



(f) Unsharp masking - filter size 4" 


Fig. B.5: Detecting shell S17 by various methods. MegaCam image. 
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